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Abstract. 

LISA is a joint space mission of the ESA and NASA for detecting low frequency 
gravitational radiation in the band 10 — 5 — 1 Hz. In order to attain the requisite sensitivity for 
LISA, the laser frequency noise must be suppressed below the other secondary noises such as 
the optical path noise, acceleration noise etc. This is achieved because of the redundancy in the 
data, more specifically, by combining six appropriately time-delayed data streams containing 
fractional Doppler shifts - time delay interferometry (TDI). The orbits of the spacecraft are 
computed in the gravitational field of the Sun and Earth in the Newtonian framework, while 
the optical links are treated fully general relativistically and thus, effects such as the Sagnac, 
Shapiro delay, etc. are automatically incorporated. We show that in the model of LISA that 
we consider here, there are symmetries inherent in the physics, which may be used effectively 
to suppress the residual laser frequency noise and simplify the algebraic approach to TDI. 



PACS numbers: 95.55.Ym, 04.80. Nn, 07.60.Ly 



Time- delay interferometry and the relativistic treatment of LISA optical links 



2 



1. Introduction 

LISA - Laser Interferometric Space Antenna - is a proposed mission which will use coherent laser 
beams exchanged between three identical spacecraft forming a giant (almost) equilateral triangle 
of side 5 x 10 6 kilometres to observe and detect low frequency cosmic GW [TJ. 

In ground based detectors the arms are as symmetrical as possible so that the laser light 
experiences nearly identical delay in each arm of the interferometer which reduces the laser 
frequency /phase noise at the photodetector. However, in LISA, the lack of symmetry will be 
much larger than in terrestrial instruments. Laser frequency noise dominates the other secondary 
noises, such as optical path noise, acceleration noise by 7 or 8 orders of magnitude, and must 
be removed if LISA is to achieve the required sensitivity of h ~ 10~ 22 , where h is the metric 
perturbation caused by a gravitational wave. In LISA, six data streams arise from the exchange 
of laser beams between the three spacecraft approximately 5 million km apart. These six streams 
produce redundancy in the data which can be used to suppress the laser frequency noise by the 
technique called time-delay interferometry (TDI) in which the six data streams are combined 
with appropriate time-delays [2]- 

This work was put on a sound mathematical footing by showing that the data combinations 
constituted an algebraic structure; the data combinations cancelling laser frequency noise formed 
the module of syzygies over the polynomial ring of time-delay operators PJ. The module was 
obtained - that is its generators were obtained - for the simple case of stationary LISA in flat 
spacetime. These were the so-called first generation TDI. However, LISA spacecraft execute a 
rotational motion, the arm-lengths change with time and the background spacetime is curved, all 
of which affect the optical links and the time-delays. The rotation gives rise to the Sagnac effect 
which implies that the up-down optical links are unequal, the arm-lengths or the time-delays 
change with time - flexing of arms. These effects cannot be ignored if the laser frequency noise 
is to be effectively cancelled. 

In this paper, we compute the orbits of spacecraft in the Newtonian framework where the 
Earth's gravitational field is also taken into account. The base orbits we take to be Keplerian 
in the gravitational field of the Sun only. On these base orbits, we linearly superpose the 
perturbative effect of the Earth's gravitational field. We choose the Earth over Jupiter because 
(i) the Earth perturbs the Keplerian orbit in resonance, resulting in a secular growth of the 
perturbations and, (ii) Jupiter's effect is less than 10% of that of the Earth's on the flexing 
and hence not dominant. The perturbative analysis is carried out within the Clohessy- Wiltshire 
(CW) [4] framework. Further, an extension of the previous algebraic approach is proposed for 
the general problem in which the time-delay operators in general do not commute; this leads 
to the second generation TDI and imperfect cancellation of laser frequency noise. However, we 
show that there are symmetries in the physical model which can simplify to some extent the 
totally non-commutative problem. These computations will be useful in the development of a 
LISA simulators, the LISACode for instance [5]. 

2. The spacecraft orbits in the Sun's and Earth's field 

The Keplerian orbits, the orbital motion in the gravitational field of the Sun only are chosen so 
that the peak to peak variation in armlengths is the least ~ 48000 km, see [6]. We summarise the 
results of paper [$] below. We choose the Sun as the origin with Cartesian coordinates {A, Y, Z} 
as follows: The ecliptic plane is the X — Y plane and we consider a circular reference orbit of 
radius R — 1 A. U. centred at the Sun. Let S — 5a/8 where a = L n /2R and L ~ 5, 000, 000 km 
is a constant representing the nominal distance between two spacecraft of the LISA configuration. 
We choose the tilt of the plane of the LISA triangle to be S = it/ 3 + do which has been shown to 
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yield minimum flexing of the arms. We choose spacecraft 1 to be at its lowest point (maximum 
negative Z) at t = 0. This means that at this point, Y = and X ~ R(l — e). The orbit of 
the first spacecraft is an ellipse with inclination angle eo, eccentricity e and satisfying the above 
initial condition. 

From the geometry, eo and e are obtained as functions of 5, 

a sin S 

tan €q 



a cos 5 + sin(7r/3) 

1/2 

-1. (1) 



1 4 2 4 
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The equations for the orbit of spacecraft 1 are given by: 
Xi = i?(cos tpi — e) cos eo, 



Y\ = R\/l — e 2 sin -01, 

Zi = — _R(cos0i — e) sineo- (2) 

The eccentric anomaly ipi is implicitly given in terms of t by, 

ipi — e sin 0i = fit — 0o , (3) 

where t is the time and ft is the average angular velocity and 0o the initial phase. The orbits 
of the spacecraft 2 and 3 are obtained by rotating the orbit of spacecraft 1 by 27r/3 and 4-7r/3 
about the Z— axis; the phases ip2, i>3, however, must be adjusted so that the spacecraft are at a 
distance ~ Lq from each other. The orbital equations of spacecraft k — 2,3 are: 

X k = Xi cos a k - Yi sin a k , 

Yfe = X\ sin (T fc + Yi cos a k , 

Z k =Z x , (4) 

where a k = (k — 1) with the caveat that the ipi is replaced by the phases ipk, where they are 
implicitly given by, 

ipk - esin0 fe = fit - (Tfc - 0o- (5) 

These are the exact (Keplerian) expressions for the orbits of the three spacecraft in the Sun's 
field. 

The Earth's field is now included perturbatively using the CW framework. The CW frame 
is chosen as follows: We take the reference particle to be orbiting in a circle of radius R with 
constant Keplerian angular velocity Q. Then the transformation to the CW frame {x, y, z} from 
the barycentric frame {X, Y, Z} is given by, 

x = (X - R cos fit) cos fit + (Y - R sin fit) sin fit , 

y = - (X - R cos fit) sin fit + (Y - R sin fit) cos fit , 

z = Z. (6) 

The a; direction is normal and coplanar with the reference orbit, the y direction is tangential and 
comoving, and the z direction is chosen orthogonal to the orbital plane. Linearised dynamical 
equations for test-particles in the neighbourhood of the reference particle are easily obtained. 
Since the frame is noninertial, Coriolis and centrifugal forces appear in addition to the tidal 
forces. With the help of the CW formalism, it is easy to see that to the first order in a (or 
equivalently e) there exist configurations of spacecraft so that the mutual distances between 
them remain constant in time. The flexing appears only when we consider second and higher 
order terms in a. In fact in [6] we have shown that the second order terms describe the flexing 
of LISA's arms quite accurately as compared to the exact Keplerian orbits. 
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The CW equations for a test particle are given by: 

x- 2Qy- 30 2 x = 0, 
y + 2VLx = 0, 

z + n 2 z = 0. (7) 

We choose those solutions of Eq. Q which form an equilateral triangular configuration of side L 
(such solutions exist). For the fcth spacecraft we have the following position coordinates: 

x k = - ^po cos(fW - cr fe - 0o ) , 
y k = po sin(fit -a k - O ) , 

z k = —p cos(ttt -a k - <fo) , (8) 

where po = L /V3. Also at t = the initial phase of the configuration is described through <p . 
In this solution, any pair of spacecraft maintain the constant distance L between each other. 

LISA follows the Earth 20° behind. We consider the model where the centre of the Earth 
leads the origin of the CW frame by 20° - thus in our model, the 'Earth' or the centre 
of force representing the Earth, follows the circular reference orbit of radius 1 A. U. Also 
the Earth is at a fixed position vector r® = (x©,?/®,z®) in the CW frame. We find that 

x B = -R(l - cos 20°) 9 x 10 6 km, y B = R sin 20° - 5.13 x 10 7 km and z® = 0. The 

acceleration field a due to the Earth at any point r (in particular at any spacecraft) in the CW 
frame is given by: 

where M® ~ 5.97 x 10 24 kg is the mass of the Earth and G = 6.67 x 10~ n kg _1 m 3 sec~ 2 Newton's 
gravitational constant. 

In order to write the CW equations in a convenient form we first define the small parameter 
e in terms of the quantity w® = GM©/d®, where d® = |r®| is the distance of the Earth from the 
origin of the CW frame; d® ~ 5.2 x 10 7 km which is more than 50 million km. So when deriving 
the forcing term we make the aprroximation |r — r®| w d®, that is, we neglect |r| compared to 
si®. It will turn out that the flexing due to the Earth is small so that this approximation is not 
unjustified. We define e = u^/fl 2 ~ 7.16 x 1CP 5 which is the just the ratio of the tidal forces 
due to the Earth and the Sun. The CW equations including the Earth's field take the form: 

x - 2Qy - SD, 2 x + eVL 2 (x - x ffi ) = , 
y + 2rix + en 2 (y - y e ) =0, 

z + n 2 {l + e)z =0. (10) 

Note that the compounded flexing due to the combined field of Earth and Sun is a nonlinear 
problem; it is infact a three body problem. We however solve this problem approximately. 
Assuming that both effects are small we may linearly add the flexing vectors due to the Sun and 
Earth; that is, add the perturbative solutions obtained from Eqs.0 and lfT7j|) : the nonlinearities 
appear at higher orders in a and e. These would modify the flexing but we may neglect this 
effect because of the smallness. We find that the flexing produced by the Earth is of the order of 
1 or 2 m/sec upto the third year, just about 40 % of that due to the Sun. But, as shown in [6] the 
flexing produced by the Sun's octupole field is nearly exact to that produced by the Keplerian 
orbits. Thus we may do better by just adding the flexing vector produced by the Earth to the 
Keplerian orbit of the relevant spacecraft. 
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We then seek perturbative solutions to Eq. I|10p to the first order in e. We write, 
x = xq + ex\, y = yo + ey±, z — zq + ez\ where xo,yo, zq are solutions at the zeroth order given 
by Eq.|(8]). We put at = (or equivalently include it in <fo) in these solutions for simplifying the 
algebra. 

With the initial conditions: x\ = y\ = Z\ = X\ = y% = z\ = at t = 0, we have the results: 

5 

x\ = — Pe cos(f2t — + + 2y©£l£ — 2pQCOS(j)o + -poQt sin(fif — 0o) > 

3 

yi = 2 i o e [sin(ftt - <j> m ) + sin0 e ] - -p [sm(Qt - cj> ) + sin0 o ] 

5 3 
+ -p Qt cos(flt ~ 4> ) - r2i(2x e - 3p cos fa) - - fl 2 t 2 y e , (11) 

where, 

5 

P% = (x@ - 2 Po cos O ) 2 + (2y ffi - 7 Po sin O ) 2 , 
- jpo sin 0o 

tan0 e = ^ — . 12 

xq - 2p cos (pa 

The z equation can be exactly integrated and used directly to obtain the flexing. However, we 
can also expand this solution to the first order in e and the result is: 

z\ = ^j-pol^sinfiicos^o — (fltcosilt — sin fit) sin^o] . (13) 

As argued before, we add the perturbation given by eri = e(xi, yi,zi) to the Keplerian orbit of 
each spacecraft. Next we compute the optical links. 



3. The six optical links 

The time-delay that is required for the TDI operators needs to be known very accurately - at least 
to 1 part in 10 8 , that is, to about few metres - for the laser frequency noise to be suppressed. In 
order to guarantee such level of accuracy, we numerically compute the optical links or the time- 
delay. This approach is guaranteed to give the desired accuracy or even better accuracy than 
what is required. We numerically integrate the null geodesies followed by the laser ray emitted 
by one spacecraft and received by the other. This computation is performed in the barycentric 
frame, and taking into account the fact that the spacetime is curved by the Sun's mass only 
(the Earth's contribution is about 5 orders of magnitude less). The computation here is further 
complicated by the fact that the spacecraft are moving in this frame of reference and the photon 
emitted from one spacecraft must be received by the other spacecraft. We use the Runga-Kutta 
numerical scheme to integrate the differential equations describing the null geodesies. But since 
the end point of the photon trajectory is not known apriori, an iterative scheme must be devised 
for adjusting the parameters of the null geodesic, in order that the worldlines of the photon and 
the receiving spacecraft intersect. We have devised such a scheme based on the difference vector 
between the photon position vector and receiving spacecraft position vector. The six optical 
links Lij have thus been numerically computed with sufficient accuracy required for TDI. The 
code gives results accurate to better than 10 metres - most of the time better than 10 -2 metres 
- except in a window of about half an hour when the error exceeds this value and becomes 
unacceptably large. We display the results in Figures [T] and [2] for O = 0. More details may be 
found in [8]. 

Figure [T] shows all the six optical links in the combined field of the Sun and Earth. We also 
need to estimate the variation in armlength which is important for the TDI analysis to follow. 
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Figure 1. The figure shows the variation in the six optical links of the LISA model for three 
years for r/>o = 0. The lengths are given in metres. 
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Figure 2. The rate of change of armlengths for the six links is shown in units of m/sec for 
4>o = 0. This rate of change is less than 6 m/sec upto the second year and increases to a 
maximum of about 8 m/sec in the third year. 

Figure [2] shows the rate of change of the six optical links as a function of time over a period 
of three years. We find that in the optimised model of LISA configuration, this rate of change 
is less than 4 m/sec. if we just consider the Sun's field. Including the Earth's field the flexing 
still remains < 6 m/sec in the first two years and increases to < 8 m/sec in the third year. 
Earlier estimates were ~ 10 m/sec. These numerical estimates are most crucial for their effect 
on residual laser frequency noise in the TDI. 

4. Time-delay interferometry 





In order to cancel the laser frequency noise, time-delayed data streams are added together in 
which an appropriate set of time-delays are chosen. In general the time-delays are multiples 
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of the photon transit time between pairs of spacecraft. In [3] a scheme based on modules over 
commutative rings was given where the module of data combinations cancelling the laser noise was 
constructed. This fully cancels the laser frequency noise for stationary LISA. There are only three 
delay operators corresponding to the three armlengths and the time-delay operators commute. 
This scheme can be straight forwardly extended to moving LISA [Zj, where, now because of 
Sagnac effect, the up and down optical links have different armlengths but the armlengths are 
still constant in time. In this case, there are six delay operators corresponding to the six optical 
links and they commute. These are the modified but still first generation TDI. However, for 
LISA the armlengths do change with time - flexing of the arms - and the first generation TDI 
modified or otherwise do not cancel of the laser frequency noise sufficiently. 

4-1. Equations for the polynomial vector cancelling laser frequency noise 

We follow the notation and conventions of [7] and [3] which are the simplest for our purpose. 
The six links are denoted by U % , V 1 , i = 1, 2, 3. The time-delay operator for the link U 2 from S/C 
1 to S/C 2 or 1 — > 2 is denoted by x in [7] and so on in a cyclic fashion. The delay operators 
in the other sense are denoted by I, m, n; the link — V 1 from 2 — > 1 by I and similarly the links 
V 2 , V 3 are defined through cyclic permutation. 

Let Ci(t) — A^j(i)/^ represent the laser frequency noise in S/C i. Let j be the delay 
operator corresponding to the variable armlength Lj(t), i.e. jCi{t) — d{t — Lj{t)). Then we 
have, 

U 1 = C x - zC 3 , 

V 1 = IC 2 - C x . (14) 

The other links in terms of Ci(t) are obtained by cyclic permutations. Also in the U l ,V l we 
have not included contributions from the secondary noises, gravitational wave signal etc. since 
here our aim is to deal with laser frequency noise only. Any observable X is written as: 

X = p i V i + q i U i , (15) 

where pi,qi,i = 1,2,3 are polynomials in the variables x,y, z,l,m,n. Thus X is specified by 
giving the six tuple polynomial vector ( Pi , qi ). Writing out the (Vi,Ui) in terms of the laser 
noises Ci(t), and in order that the laser frequency noise cancel for arbitrary Ci(t), the polynomials 
(Pi,qi) must satisfy the equations: 



Pi 


- ?l - 


V q2x 


-p 3 n 


= 0, 


Pi 






-pil 


= 0, 


p?, 


- <73 - 


\~ qiz 


- p 2 m 


= 0. 



The solutions to these equations as realised in earlier works are important, because they consist 
of polynomials with lowest possible degrees and thus are simple. Since these are linear equations 
they define a homomorphism of modules and the solutions themselves form a module - the 
module of syzygies over the polynomial ring Q{x,y,z,l,m,n), where Q is the field of rational 
numbers and x, y, z, I, to, n play the role of indeterminates. In general, the variables (operators) 
x,y, z,l,m,n do not commute and hence the order of the variables is important. However, if 
we assume in a simple model that the arms do not flex, then the operators commute, and the 
generators of the module have been found via Grobner basis methods [31 E|- 

However, when the arms flex, the operators no longer commute. If we operate on C(t) with 
operators j and k in different orders, it is easily seen that jk ^ kj. A combinatorial approach 
has been adopted in [9] to deal with the totally non-commutative case. However, our aim here 
is to estimate the level of the non-commutativity of these operators in the context of our LISA 
model and use the symmetries to simplify the algebraic approach. 
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4-2. The algebraic approach and symmetries 

The level of non-commutativity can be found by computing commutators which occur in several 
of the well known TDI observables like the Michelson, Sagnac etc. We find that given our 
model of LISA, we require to go only upto the first order in L; we find for our model L ~ 10~ 6 
metres/sec 2 and thus even if one considers say 6 successive optical paths, that is, about At ~ 100 
seconds of light travel time, At 2 L ~ 10~ 2 metres. This is well below few metres and thus can 
be neglected in the residual laser noise computation. Moreover, 1? terms (and higher order) can 
be dropped since they are of the order of < 10~ 15 (they come with a factor 1/c 2 ) which is much 
smaller than 1 part in 10 s . The calculations which follow neglect these terms. 

Applying the operators twice in succession and dropping higher order terms as explained 
above, 

k 2 kiC = C(t - L kl (t - L k2 ) - L fca ) , 

« C{t - L ki - L k2 ) + L k2 L kl C(t - L kl -L k2 ). (17) 

The above formula can be easily generalised by induction to n operators. 

We now turn to the commutators of the operators. The term in C cancels out; only the C 
term remains. We list below a few of the commutators: 

jk - kj = LjL k - L k Lj , 

Imjk - jklm = (Li + L m )(Lj + L k ) - (Lj + L k )(Li + L m ) , 

Imnxyz — xyzlmn = (Li + L rn + L n )(L x + L y + L z ) (18) 
— (L X + Ly + L z )(L l +L m + L n ). (19) 
We observe the following approximate symmetries in our model: 

L x ~ L[, L y f=a L m , L z f=a L n , (20) 
which also implies (this combination occurs in the Sagnac observables), 

L x + L y + L z fa Li + L m + L n . (21) 

Infact in our model, \(L X + L y + L z ) — (L; + L m + L n )\ < 1 m/sec and \L X — Li\ < 0.8 m/sec 
upto the first three years in our model. The same is essentially true for the pairs of links y, m 
and z, n. Thus these pairs of operators essentially commute. Thus, we are not dealing with a set 
of totally non-commuting variables, but with an intermediate case. 

In addition to these approximate symmetries there are other commutators which vanish 
'identically' (after dropping terms in L 2 and L and higher order). It can be easily verified that 
commutators of the form, [x\X2---x n , yiy2---yn], n>2 vanish identically, when the y\,yi, ■•■iVn 
are permutation of the operators x\, x n , and the x k for a given k represents one of the 

delay operators x, y, z, I, to, n. Thus the original non-commutative ring Q(x, y, z, I, to, n) can be 
quotiented by the ideal U which is generated by the vanishing and approximately vanishing 
commutators: U = {[x, I], [y, to], [z, n]; [xy, yx], ...}. The quotient ring Q/U is clearly much 
smaller and simpler and the solution to Eq. (fT6]) is sought for polynomials in this quotient 
ring. The solution set of polynomial vectors {pi,qi) still form a module over Q/U. The future 
goal is to 'construct' this module. 

4-3. Residual laser frequency noise in the Sagnac observable and symmetries 

By the time LISA flies the expectations are for the laser frequency noise estimate to reduce to say 
Av ~ 10Hz/\/Hz. If we divide this number by the laser frequency v o ~ 3 x 10 14 Hz, we obtain 
the noise estimate in the fractional Doppler shift C with the power spectral density (PSD): 

Sc(f) = (\C(f)\ 2 ) - 10- 27 Hz" 1 , (22) 
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where C(/) is the Fourier transform of C(t). Then by differentiating C, the PSD of the random 
variable C is just S d (f) = 4tt 2 / 2 5 c (/) Hz. 

The modified Sagnac first generation TDI observable a is given by the polynomial vector in 
the form (pi,qi) by: 

a = (k, kI, Klm,r),r)zy,riz) , (23) 

where k = 1 — zyx and 77 = 1 — Imn. If the variables x,y, z,l,m,n commute then the laser 
frequency noise is fully cancelled. However, if they do not commute, there is a residual term. It 
can be computed as: 

AC = a x Cx + a 2 C 2 + a 3 C 3 ■ (24) 

We find that a 2 = a 3 = and ct\ = [zyx, Imn] and so by Eq. lfl9|) : 

At(i) = ^ + L y + L»)(L| + L m + L„) - (L, + L m + L n )(L x + L y + L z )}, (25) 

and thus AC = AtC\. Because the vary during the course of an year the At also varies 
during the year and so also the amplitude of the random variable AC. Thus the PSD of AC is: 

S AC (f;t) = 4TT 2 At(t) 2 f 2 S c (f). (26) 

This noise must be compared with the secondary noise pQ. However, because we are considering 
the modified TDI Eq. lf23|) . there are extra factors k and r\ which do not appear in the 
corresponding first generation TDI. These factors introduce an additional multiplicative factor, 
namely, 4 sin 2 (37r/L ) in the secondary noise PSD which leaves the SNR unchanged but must be 
considered when it is compared with the residual laser frequency noise given in Eq. I|26p . Thus, 

S a (f) = 4sin 2 (37r/Lo){[8sin 2 3^/L + 16sin 2 7r/L ]5 acc + 65 opt } (27) 

where S acc = 2.5 x 10- 48 (//lHz)- 2 Hz" 1 and S opt = 1.8 x 10- 37 (//lHz) 2 Hz" 1 . In the Figured 
we plot S a (f) and Sac(/;*) at three epochs an year apart. 

Comparison of noise PSDs 




1e-65 1 1 1 — 1 1 — 1 1 — 1 1 — 1 

1e-04 0.001 0.01 0.1 1 

fin Hz 

Figure 3. The 'top' curve shows the PSD S a (f) of the secondary noises. The straight lines 
are the PSDs of the residual noise at three epochs chosen an year apart. Clearly the residual 
laser noise is adequately below the secondary noises. 
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We see that, clearly the residual laser frequency noise is few orders of magnitude below the 
secondary noises. Since the other Sagnac variables j3 and 7 are obtained by cyclic permutations 
of the spacecraft, the residual laser noise is similarly suppressed in them. The basic reason for 
this remarkable cancellation is the symmetry inherent in the physics. Note that a\ £ U and 
hence a is an element of the module we are seeking. 

5. Concluding remarks 

We have computed in the Newtonian framework the spacecraft orbits in the combined field of 
the Sun and Earth and from this deduced the flexing of the arms of LISA by choosing the model 
which gave minimum flexing when only the Sun's field was taken into account. Now the flexing is 
no more periodic as was the case when only the Sun's field was considered. We have ignored the 
effect of Jupiter because we believe this effect to be not so dominant as that of the Earth. Writing 
the tidal parameter for Jupiter, e.j = GMj/dj, similar to e of Earth, where, Mj w 2 x 10 27 kg 
is the mass of Jupiter and dj, the distance from LISA to Jupiter, which we take on the average 
to be ~ 5 A. U., we find ej/e <~ 0.09. Moreover, Jupiter has its own periodicity pertaining to its 
orbit and therefore will not be in resonance as was the case with the Earth, and thus there will 
be no secular effect. Thus we do not expect the effect of Jupiter on flexing to dominate. Note 
that these results are valid so long as we can neglect the nonlinearities arising from higher order 
terms in e and a. 

We have computed the residual laser frequency noise in one of the important TDI variables, 
namely, the Sagnac. The residual noise is satisfactorily suppressed because of the symmetry. 
In other variables such as the Michelson this is not true and higher degree polynomials will be 
required. The algebraic approach outlined above seems promising. 

Our model of LISA is optimal (minimal flexing of arms) only in the Sun's field. Clearly this 
opens up the question of seeking an optimal model for the LISA configuration in the field of the 
Sun, Earth, Jupiter and other planets which will minimise the flexing of the arms and therefore 
the residual laser frequency noise in the modified first generation TDI. We finally remark that 
our computations here may be useful in the development of a LISA simulator. 
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